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purpose.  An  alternative  approximate  method  has  also  been 
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Summary 

To  make  best  use  of  ionospher Ic  ref lect Ions  for  very  long  range  transmissions 
of  radio  waves  in  the  HF  part  of  the  spectrum,  it  is  necessary  to  radiate 
vertically  polarized  waves  having  substantial  field  strength  in  the  directions 
near  the  horizon.  However,  as  is  well  known,  the  imperfect  conductivity  of 
typical  earth  tends  to  reduce  drastically  the  available  signal  near  the  horizon. 
The  possibility  of  increasing  the  radiation  at  such  low  angles  by  the  inter- 
action of  a horizontal  ground  wire  system  is  known  and  has  been  discussed  but 
heretofore  there  has  been  very  little  quantitative  information  available  on 
this  subject  and  no  systematic  procedures  for  the  design  of  such  ground  systems; 
in  view  of  the  substantial  cost  of  such  ground  systems,  much  better  information 
on  their  performance  is  needed. 

The  work  described  in  this  interim  report  has  as  its  basic  purpose  the 
systematic  determination  of  the  effects  that  horizontal  ground  wire  systems 
may  have  on  the  low  angle  radiation,  and,  eventually,  the  development  of 
ra^onal  procedures  for  the  design  of  such  ground  systems  for  application  in 
specific  tactical  systems.  In  the  course  of  this  study,  it  has  been  demonstrated 
that  useful  results  can  be  obtained  with  the  aid  of  an  extensive  digital  computer 
program  incorporating  the  best  available  theoretical  representat ion  (Sommerfeld 
integrals)  in  a "moment  method"  for  the  determination  of  antenna  currents  and 
fields.  It  has  also  been  shown  (by  comparison)  that  useful  results  may  be 
obtained  from  a relatively  simple  digital  computer  program  that,  unfortunately, 
must  be  based  on  assumed  current  distributions  in  the  ground  wires. 

The  amount  of  low  angle  radiation  enhancement  that  is  feasible  by  design 
of  such  ground  wire  systems  depends  on  earth  conductivity  and  frequency.  It 
ranges  from  just  a fewpercent  at  the  lower  end  of  the  HF  band  to  25%  or 
more  at  the  high  end  of  the  HF  band  in  a region  having  poor  earth  conductivity. 


ill 


1.  Introduction 


For  antennas operat  ing  In  the  HF  range  or  lower,  antenna  performance  suffers 
unless  the  typical  poor  earth  conductivity  found  at  most  sites  is  enhanced  by 
a system  of  conducting  wires  near  or  below  the  surface  of  the  earth.  This 
ground  wire  system  is  usually  essential  for  impedance  (l.e.,  lowered  input 
resistance  losses)  control,  but  it  might  also  be  helpful  in  radiation  pattern 
control,  especially  to  alleviate  the  well  known  problem  of  the  null  at  the 
horizon  found  w I th  vertically  polarized  antennas.  This  report  deals  with  only 
the  latter  problem;  it  describes  investigations  into  the  extent  to  which  a ground 
wire  system  may  be  used  to  enhance  the  radiation  pattern  of  some  simple  vertically 
polarized  arrays. 

Evaluation  of  the  effects  of  a ground  screen  on  or  in  an  actual  earth 
by  standard  boundary  value  solutions  of  the  field  equations  is  impractical 
at  the  present  time.  However,  over  the  past  10  to  12  years,  a method  of 
calculation  of  current  distribution.  Impedance,  radiation  and  scattering 
patterns  of  wire  like  structures  has  been  highly  developed.  A number  of 
computer  programs,  based  on  the  method  (that  has  come  to  be  known  as  the 
method  of  moments)  have  been  developed  to  handle  wire-1  Ike  structures  of 
considerable  complexity  and  variety.  In  the  main,  the  programs  treat  the 
antennas  acting  only  In  free  space,  or  at  most,  in  the  vicinity  of  a perfect- 
ly conducting  ground.  A group  at  Lawrence  Livermore  Laboratories,  however, 
modified  and  extended  a program'of  this  type  (WAMP)  to  cover  wire  structures 
in  or  near  a plane  homogeneous  earth  having  arbitrary  conductivity.  They 
called  the  program  WF-LLL2A  and  described  it  in  the  "Electromagnetic  Computer 
Code  Newsletter"  Volume  2 No.  1 dated  April  1,  1975  and  in  UCRL-51821  ^ 
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"Fortran  Subrountlnes  for  the  Numerical  Integration  of  Sommerfeld  Integrals 
Unter  Anderem",  by  Lager  and  Lytle,  dated  May  21,  1975-  A group  at  Purdue 
had  been  pursuing  a similar  development  and  upon  learning  of  the  Livermore 
work  decided  to  try  to  adopt  and  adapt  the  WF-LLL2A  program  for  the  purpose 
of  doing  ground  screen  calculations,  although  neither  the  WF-LLL2A  program 
nor  the  modifications  of  It  had  yet  been  tested  or  evaluated  for  the  purpose. 

But  the  basic  idea  was  to  restrict  potential  ground  screen  systems  to 
wire  structures,  and  to  regard  the  wires  in  the  ground  screen  as  part  of 
the  antenna  in  or  near  the  earth.  By  and  large,  this  meant  that  the  fields 
due  to  the  elemental  currents  flowing  in  the  antenna  wires  had  to  be  evaluated 
with  the  Sommerfeld  integrals,  or  suitable  approximations  to  them.  It  turned 
out  that  none  of  the  existing  codes  would  converge  properly  when  the  wire 
elements  were  close  to  the  earth  and  to  each  other.  Locating  the  source  of 
the  difficulty  and  correcting  it  involved  considerably  more  time  and  effort 
than  had  been  anticipated. 

2.  Background  of  the  Wire  Antenna  Program 

The  computer  programs  assume  the  structure  consists  of  thin  wires  and 
essentially  solve  an  electric  field  integral  equation,  numerically,  for  the 
currents  that  result  with  a stated  applied  eiectric  field  (or  voltage).  The 
program  WF-LLL2A  employs  a three  term  current  expansion  function  for  the  cur- 
rent in  each  subsection  of  wires. 

The  essence  of  the  method  is  to  replace  an  integral  equation  for  the 
eiectric  field  by  a set  of  linear  equations  relating  the  electric  field  at 
a set  of  discrete  points  to  a set  of  current  levels  that  produce  the  fields 


at  the  set  of  points.  The  coefficients  in  the  set  of  linear  equations  have 
the  character  of  impedance  coefficients  and  the  matrix  relating  the  electric 
fields  and  current  level s is  often  called  the  impedance  matrix,  or  system  matrix 
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or  system  Impedance  matrix.  The  task  here  Is  two  fold  1)  first  compute  the 
elements  of  the  system  impedance  matrix,  which  means  calculate  the  field  due 
to  a unit  amplitude  current  segment  in  the  appropriate  environment  at  each  of 
the  field  points  2)  having  obtained  the  Impedance  matrix,  solve  the  resulting 
, set  of  linear  equations  for  the  currents  (l.e.,  Invert  the  matrix). 

The  difficulty  of  applying  this  method  to  the  calculation  of  the  effects 
of  wire  ground  screens  on  antenna  arrays  is  two-fold.  First,  it  is  difficult 

i 

to  calculate  the  field  produced  by  current  elements  at  nearby  points  when  the 

•»% 

current  element  is  near  the  surface  of  an  imperfectly  conducting  earth.  Second, 
with  the  ground  wire  system  plus  the  antenna  wires  the  number  of  segments 
(i.e.,  the  number  of  equations  in  the  linear  set)  into  which  the  wire  con- 
figurations must  be  broken  is  very  large.  Thus,  the  basic  problem  is  one  of 
time  and  storage  space  in  the  digital  computer. 

The  WF-LLL2A  has  available  within  it  four  methods  for  the  calculation  of 
the  fields  in  the  presence  of  ground.  Two  of  these  are  based  on  plane  wave 
reflection  coefficients  to  account  for  the  reflected  fields.  The  fields  that 
would  result  from  a perfectly  conducting  ground  (i.e.,  an  Image ) are  modified 
according  to  the  actual  conductivity  and  dielectric  constant  of  the  ground. 

A third  method  makes  use  of  the  well  known  formulas  developed  by  Norton; 
this  method  is  used  when  the  distances  and  ground  parameters  are  in  the  proper 
range  for  the  Norton  formulas.  The  fourth  method  Incorporated  is  the 
numerical  evaluation  of  the  Sommerfeld  integrals,  and  for  this  there  are  two 
options  available,  the  Hankel  function  form  and  the  Bessel  function  form,  with 
the  choice  between  the  two  being  made  on  the  basis  of  position  and  parameter 
values.  The  Integrals  are  evaluated  by  performing  a contour  integration  In 
the  complex  plane,  numerically,  using  an  adaptive  Romberg  Integration. 


4 


3.  Modifications  for  Calculation  of  Effects  of  Ground  Wire  Systems 

In  the  original  program,  for  points  nearby  thin  wires  located  close  to  the 
surface  of  the  earth,  the  Sommerfeld  integral  calculations  referred  to  above 
are  extremely  time  consuming,  if  indeed  they  converge  properly  at  all. 

It  was  recognized  from  the  outset  that  computer  time  would  be  a problem, 
so  the  program  was  reorganized  in  several  respects.  First,  since  the  most 
time  consuming  calculations  are  nearby  "self"  and  "mutual"  impedance  coefficients, 
the  matrix  was  partitioned  so  that  these  time  consuming  calculations,  which 
it  was  observed,  actually  need  to  be  done  only  once  for  a given  size  wire, 
frequency,  and  ground  parameters  set,  could  be  calculated  once  and  for  all 
and  fed  in  thereafter,  when  successive  runs  were  being  made  for  the  purpose 
of  evaluation  of  the  effects  of  different  sizes  and  shapes  of  ground  screens. 

Other  programming  techniques  were  developed  so  as  to  conserve  the  space 
available  in  memory  by  having  in  the  fast  memory  only  those  parts  of  the 
program  needed  for  the  particular  calculation  being  done. 

A basic  difficulty  was  encountered,  however,  in  the  use  of  the  Sommerfeld 
integral  portion  of  the  program,  when  it  was  discovered  that  for  interesting 
wire  sizes,  ground  parameters  and  positions,  a time  of  1 048  seconds  in  the 
CDC  6500  computer  (one  of  the  "priority"  time  limits  set  by  The  Purdue  Computing 
Center)  was  insufficient  to  calculate  the  self  impedance  coefficients  for  even 
a very  simple  single  wire  structure.  Moreover,  the  indications  were  that 
even  a significantly  larger  amount  of  time  would  also  be  insufficient.  It 
appeared  then  that  some  of  our  mod i f lc iat I ons  to  the  calculation  would  be 
essential  and  a step-by-step  replacement  and  checking  of  various  parts  of  the 
existing  Sommerfeld  integration  packages  was  undertaken.  Suffice  it  to  say, 
that  although  most  of  the  modifications  decreased  the  time  necessary  for  the 


calculation,  they  did  not  decrease  it  enough  to  be  really  useful  and  even  so, 
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some  were  giving  up  accuracy  in  the  calculation.  A very  considerable  amount 
of  time  was  spent  in  tracking  down  the  real  source  of  the  difficulty,  and  in 
trying  to  determine  the  accuracy  of  the  modified  systems. 

To  explain  the  difficulty  involved  in  carrying  out  the  numerical 
integration,  consider  the  expression  for  the  radial  component  of  the  electric 
field  in  the  medium  above  the  ground  due  to  a horizontal  electric  dipole 
moment  p oriented  along  the  x-axis  with  notation  as  in  the  WF-LLL2A  program: 
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-Y2(h+z)  . . 

U22  = H^;(Xp)XdX  (5) 

+ y2 

2 i . 

and  k , = [id  p £ £ , - I'top  a,]2,  k„  = w/p  e are,  respectively,  the  complex 
l o o r I o l z o 0 

2 2 i 

propagation  constants  for  medium  1 (ground)  and  2 (air);  and  = [X  - kj]  , 

Y2  = [X2  - k2]^.  [To  avoid  possible  confusion,  we  have  assumed  e^1  type  of 

time  dependence.  This  change  dictates  the  use  of  the  Hankel  function  of  the 

second  kind  and  a different  choice  in  branch  cuts.  All  other  notations  are 

the  same  as  those  of  Lager  and  Lytle.] 

To  evaluate  the  elements  of  the  impedance  matrix,  the  tangential  component 

of  the  electric  field  at  the  surface  of  the  wire  is  integrated  with  respect 

to  the  dipole  moments  situated  along  the  center  of  the  wire.  For  example, 

consider  the  self  impedance  of  a segment  of  wire  of  radius  a and  length  £. 

The  field  at  x=0,  y=a,  z=h,  is  integrated  along  a path  described  by  the 

l l 

conditions  -y  < x < y,  y=0,  z=h,  as  indicated  in  Fig.  1. 

First,  focus  attention  to  the  first  part  of  the  expression  (GF I ELDS  and 
GEVALA) . When  the  wire  is  several  radii  above  ground,  h>>a,  is  a rela- 

tively smooth  function  of  x while  G22  has  a peak  at  x=0  with  a peak  value  of 

the  order  of  a When  the  wire  is  close  to  the  ground,  h~a,  both  G2j  and 

^ 2 

^22  ^ave  Pea^s  at  x=0,  and  — J (G 22  “ G i)  ^as  3 s^arP  Peak  with  a peak  value 

“5 

of  the  order  of  a . In  the  original  WF-LLL2A  program,  the  integration  is 
performed  via  an  adaptive  Romberg  method.  Due  to  the  sharp  peak,  the  integra- 
tion is  extremely  time  consuming  for  thin  wires  close  to  ground.  Actually, 
the  sharp  peak  may  be  smoothed  by  partial  integration.  We  replaced  this 
portion  of  program  (GFIELDS  and  GEVALA)  by  two  closed  form  expressions  and  an 
integration  of  an  integrand  with  a moderate  peak. 


Now  consider  the  second  part  of  equation  (1)  (SFIELDS  and  EVALUA2  (or 

EVALUA3)).  Double  integrals  are  involved  in  calculating  the  parts  of  this 

contribution  to  the  impedance  matrix,  and  the  evaluation  of  V ^ and 

(equat ions  (A)  and  (5))  is  particularly  time  consuming.  For  |h+z|  * |x|,  the 

exponential  functions  in  (A)  and  (5)  decay  rapidly  and  thus  V ^ and  ma/ 

be  evaluated  routinely.  However,  when  |h+z|  <<  | X | , the  integrands  of  (4)  and 

(5)  behave  like  |x|^  e for  | Xp  | » 1.  Thus,  the  integrands  decay 

extremely  slowly  and  the  integration  has  to  be  carried  outat  least  to  | X | >> 

p In  addition,  the  contour  has  to  be  deformed  in  the  complex  X-plane  so 

(2) 

as  to  avoid  the  poles  and  branch  cuts.  The  Hankel  function  y ( Xp ) with 
complex  argument  Xp  is  extremely  difficult  to  evaluate  particularly  for  the 
case  p <<  1 . Since  our  concern  is  mainly  for  the  region  where  p <<  1 , h+z  « 1, 
which  means  that  the  time  retardation  is  negligibly  small,  we  replace  (A)  and 
(5)  by  their  quasi-static  approximations,  i.e.,  y^  ~ ~ X.  With  such 

approximations,  (4)  and  (5)  may  be  evaluated  in  closed  forms.  Our  calculation 
shows  that  in  the  region  of  interest,  these  closed-form  expressions  are  very 
accurate.  In  our  version  of  SFIELDS,  the  quasi-static  approximation  is  used 
for  p <_  5a . 

In  these  ways  we  reduce  the  time  required  for  self  and  nearby  mutual 
impedance  terms  to  manageable  proport  ions  and  have  thus  arrived  at  our  "best 
method  to  date".  We  will  refer  to  this  modified  program  as  WF-LLL2AP. 

4.  Alternative  Method  of  Calculation 
An  approximate  but  very  fast  method  of  evaluating  the  effect  of  ground 
wires  on  the  radiation  patterns  of  vertical  antennas  over  actual  grounds  may  be 
based  on  the  reciprocity  theorem.  Fig.  2 shows  an  antenna  array  on  an  imperfect 
earth.  The  current  distribution  in  the  array  is  defined  to  be  J and  is  assumed 
known  (in  ntost  cases  a simple  sinusoidal  distribution  along  the  wires  is  a 
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sufficiently  good  approximation).  The  open  circuit  voltage  at  this  input  of  the 
array  can  be  found  from  the  integral 


Voc  (a  >♦> 


J • E dv 
-'array 


assuming  a unit  current  at  the  input  on  transmitting.  Iri  this  expression 
E is  the  electric  field  of  a plane  wave  arriving  at  angle  (6,<f>)  in  the  presence 
of  the  imperfect  earth  (and  of  course  evaluated  at  the  positions  of  the 
currents  in  the  elements  of  the  array),  but  in  the  absence  of  the  array. 

Vqc  ( 6 , d> ) is  then  the  receiving  voltage  pattern  of  the  array  over  an  imperfect 

ground  and  by  reciprocity,  it  is  also  the  transmitting  pattern.  Therefore,  we 

3* 

might  as  well  write  the  result  in  terms  of  the  distant  E field 


E(e,<>)  = 


E.  • J„  dv 
i nc  t 


where  is  the  current  distribution  (as  a transmitter)  on  the  antenna  system 
when  it  is  transmitting.  The  advantage  of  this  formulation  is  that  except  for 
grazing  angles,  the  field  E.nc  may  be  determined  with  plane  wave  reflection 
coefficients  (for  grazing  angles  the  Sommerfeld  solutions  or  equivalent  should 
be  used) . 

Moreover,  if  the  ground  system  consists  of  a set  of  ground  wires,  it  is 
both  consistent  and  practical  to  regard  the  currents  in  the  ground  wires  as 
a part  of  the  transmitting  antenna  current  distribution,  that  is 

dt  dvert  + dHGW 

where  the  two  terms  are,  respectively,  the  current  density  in  the  vertical 
elements  and  (J^q^)  the  current  density  in  the  horizontal  ground  wire.  Then 
the  radiated  field  is 


E(0,<t>)  = E . • J dv  + 

J inc  vert 


E. 


i nc 


JHGW  dv 


(Weeks,  W.  L.,  Antenna  Engineering,  pg.  2b  equation  23) 
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It  Is  to  be  expected  that  Jvert  will  not  change  significantly  in  distribution 
whether  there  is  a ground  wire  system  or  not,  so  with  unit  current  at  the  input 
In  both  cases: 

e (e , <*►)  ] 

no  ground  screen 

Thus,  it  is  possible  to  define  an  enhancement  of  the  radiated  field  due 
to  the  presence  of  the  ground  wire  system  as 


i nc 


vert 


dv 


E(e,<t») 

rre: 


= 1 + 


i nc 


JHGW  dv 


i nc 


J , dv 
vert 


no  ground  screen 

For  example,  if  the  ground  is  perfect  and  the  ground  system  wires  lie  on 
the  perfect  ground,  then  Ejnc 


dHGW  = ® and  no  er|hancement  '5  possible.  In 


this  case 


E(0,$)  = J E.  • J . dv 
y J inc  vert 

is  the  pattern  of  the  transmitting  antenna  array  over  a perfect  ground. 

If  we  specify  that  the  elements  are  along  the  vertical  direction  z and 

the  plane  wave  is  incident  at  angle  0 in  the  x-z  plane  with  H perpendicular  to 

the  plane  of  Incidence,  the  E,  field  has  a form  as  follows  (in  the  air  above 

Inc 

the  ground) 


E . 
z,  i nc 


ex 

— — e 

“ E0  0 


+ pH1  (e ) e 


-jBM 


— — H e 

U)  G 0 

o 


JM  jB  z 
[e 


pH1(e)  e z ] 


-JM. 


HI , 


where  0 = 0,,sin0  and  0 ■*  0^cos0  and  3^  = «)/unen;p  (0) 

x 0 z 0 0 u u 

at  angle  0.  In  the  ground  (earth) 


ref 
* i nc 


with  wave  incident 


~7~~~  (1  + pHl(0))e  YzeZ  + YxeX 
jme 


where  z Is  the  depth  to  be  Inserted  as  a negative  number,  and 
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The  quantities  A and  B in  the  expression  lHGW(L)  are  determined 
primarily  by  the  terminations.  However  if  the  attenuation  is  large,  only 
one  of  the  constants  is  important  (B,for  wires  extending  toward  positive  x) . 

In  relating  the  current  in  the  horizontal  ground  wires  to  the  current  at 
the  input  (or  reference  point)  of  the  vertical  antenna,  a difficult,  ticklish, 
but  very  practical  consideration  arises.  For  the  field  enhancement  due  to  the 
ground  wires  certainly  depends  on  the  relative  level  of  the  currents  in  the 
ground  wires;  yet  the  level  relative  to  the  current  in  the  verticals  depends 
on  the  minute  details  of  the  method  of  attachment  of  the  transmitting  structure 
to  the  ground  and  to  the  ground  wires  at  the  point  of  feed.  The  significance 
of  this  is  suggested  by  the  sequence  of  sketches  in  Fig.  3.  Sketch  a suggests 
a limiting  situation  in  which  no  current  at  all  would  flow  in  the  ground 
wire.  The  transmitter  in  £ is  connected  both  to  a perfect  ground  and  to 
a slightly  elevated  ground  wire,  in  this  case,  essentially  all  of  the  current 
would  flow  into  the  perfect  ground.  On  the  other  hand,  another  extreme  is 
indicated  in  sketch  b^  in  which  there  is  no  direct  connection  to  the  ground. 
Although  it  might  seem  at  first  sight  that  all  of  the  current  would  flow 
into  the  slightly  elevated  ground  wire,  this  would  not  in  fact  happen  for 
there  would  be  some  hard  to  define  capacitive  susceptance  to  ground  which 
would  make  the  ground  wire  current  somewhat  less  than  the  vertical  input 
current,  as  suggested  in  sketch  c_.  Sketch  d_  shows  something  more  like  an 
actual  situation  with  a ground  connection  to  the  imperfect  earth  as  well  as 
to  the  horizontal  ground  wires.  In  this  case,  the  impedance  to  earth  will 
be  significant,  complex,  and  sensitive  to  detail;  and  the  input  current  at 
the  vertical  will  thus  d?  ide  in  some  (hard  to  determine)  fashion  between 
the  horizontal  ground  wires  and  the  grounding  connection.  From  the  point 
of  view  of  trying  to  maximize  the  enhancement  of  fields  by  the  ground  wires, 
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FIG.  ? MODELS  ILLUSTRATING  THE  FACTORS  THAT  INFLUENCE 
THE  RELATIVE  LEVEL  OF  CURRENT  IN  THE  GROUND  WIRES 
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It  Is  necessary  that  the  ground  wires  have  all  of  the.  current.  This  Implies 
the  highest  possible  grounding  Impedance.  Yet  other  considerations  (safety 
and  lightning  for  example)  require  that  the  Impedance  to  ground  be  small. 

This  suggests  that  a technique  for  grounding  which  has  a high  Impedance  at 
HF  but  a low  Impedance  at  lower  frequencies  should  be  a worthwhile  step  to 
take  in  a system  design  in  which  the  ground  wires  are  being  called  upon  to 
enhance  the  radiation  pattern. 

Another  practical  consideration  that  influences  the  enhancement  possible 
by  a ground  wire  system  of  a given  size  is  the  manner  of  termination  of  the 
wires(real ly,  the  terminating  Impedance  of  the  ground  wires  at  the  ends 
remote  from  the  transmitter).  This  termination  impedance  may  be  low  (if 
connected  to  a low  impedance  ground  electrode)  or  high  (if  the  ends  are 
capped  with  insulating  material)  or  even  "matched",  if  enough  is  known  about 
the  environment  of  the  ground  wires.  Depending  on  the  actual  terminating 
Impedance,  of  course  the  "input  impedance"  to  the  ground  wire  system  viewed 
from  the  transmitter  is  different,  and  therefore  a different  relative  level 
of  current  actually  flows  Into  the  ground  wires.  In  the  case  of  earth  having 
low  conductivities,  certain  lengths  would  exhibit  resonance  type  effects 
which  would  make  the  effective  Input  impedance  to  the  ground  wires  quite 
high  (for  example,  half-wavelength  wires  left  open  ended,  or  quarter  wavelength 
wires  having  low  impedance  terminations) 

All  of  these  effects  make  It  difficult  to  determine  what  the  levels  of 
current  in  the  ground  wires  are  relative  to  those  in  the  vertical  wires. 

Thus,  except  for  displaying  the  sensitivity  of  the  low  angle  radiation 
enhancement  to  such  effects,  for  the  most  part  in  the  approximate  calculation, 
we  will  assume  a traveling  wave  distribution  on  the  ground  wire,  with  equal 
current  levels  flowing  Into  the  verticals  and  to  the  horizontal  ground  wire  system. 


■ 
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The  output  of  the  computer  program  WF-LLL2AP  can  give  guidance  as  to 
which  assumptions  or  conditions  are  likely  to  be  the  better. 

5.  Results  of  the  Calculations 

The  overall  effects  that  radiating  ground  wires  have  on  vertical  radiation 
is  suggested  by  the  radiation  patterns  displayed  in  Figs.  4 and  5.  These 
patterns  were  actually  obtained  by  the  approximate  method  described  above,  but 
Figs.  6 and  7 show  similar  results  obtained  by  the  modified  WF-LLL2AP  program 
and  the  results  are  then  expected  to  be  typical.  The  patterns  show  EQ(s) 
in  the  forward  direction  of  the  vertical  plane.  In  each  case,  the  lower 
graph  shows  the  pattern  of  the  vertical  array  with  no  ground  wires  at  all. 

The  figures  display  the  tendency  of  the  field  to  have  a null  on  the  horizon 
(0=90°)  as  expected;  Fig.  A gives  the  result  for  earth  conductivity  10 
milllmhos/m.  Fig.  5 gives  those  for  1 mi  11  imho/m.  The  upper  curves  in  each 
case  give  the  pattern  with  a ground  wire  system  consisting  of  four  horizontal 
ground  wires  at  the  surface  of  the  earth  in  the  forward  direction  only, 
extending  to  a length  of  6 m (at  10  MHz).  (The  specific  results  depend  on 
the  ground  wire  length,  and  as  will  be  seen  below;  longer  ground  wires  do 
not  necessarily  provide  more  low  angle  field  enhancement).  Both  of  the 
figures  show  that  ground  wires  enhance  the  radiation  at  all  angles,  assuming 
that  the  current  in  the  ground  wires  near  the  input  is  equal  to  that  in  the 
vertical  wires,  and  that  the  current  distribution  in  the  ground  wires  is 
that  of  a pure  traveling  wave  at  a propagation  constant  as  calculated  by 
Coleman's  formula.  It  Is  also  seen  that  the  enhancement  is  greater  with  the 
lower  earth  conductivity,  as  might  be  expected. 

The  radiation  patterns  are  all  changing  so  rapidly  at  low  angles  (6 
near  90°)  that  the  positive  effects  of  the  ground  wires,  if  any,  are  hard  to 
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judge  from  the  patterns.  For  this  reason,  the  balance  of  data  presented  does 
not  give  the  whole  radiation  pattern  but  only  the  enhancement  ratio  (field 
magnitude  with  ground  wires/field  magnitude  with  no  ground  wires)  for  some 
particular  low  angle.  For  other  low  angles,  the  enhancement  is  not  greatly 
different  although  it  tends  to  be  slightly  smaller  for  the  lower  angles. 

The  next  Figure,  8,  shows  enhancement  data  at  10  MHz,  0=75°,  for 
three  different  earth  conductivities,  as  a function  of  the  length  of  the 
horizontal  ground  wires  ( same  assumpt ions  as  above).  These  data  show  that 
for  conductivity  of  10  millimhos/m,  the  enhancement  is  about  as  high  as  it 
ever  will  be  when  the  ground  wires  are  k to  5 meters  long,  at  which  time,  the 
maximum  enhancement  possible  would  be  at  10%.  And  in  fact,  the  field 
enhancement  goes  down  slightly  for  lengths  twice  as  long  as  the  lengths  that 
give  the  maximum  effect.  If  the  conductivity  Is  5 millimhos/m,  an  enhancement 
of  about  15?  is  possible  for  ground  wires  of  length  about  6 m and  with  longer 
wires,  the  possible  enhancement  is  less,  falling  to  about  10%  at  a length  of 
about  15  meters.  For  the  still  lower  conductivity,  namely  1 mi  11  imho/m, 
more  enhancement  is  possible,  about  20%  with  a length  of  6 meters.  However, 
in  the  data  for  1 ml  11  imho/m,  a cyclic  variation  of  enhancement  level  with 
length  is  apparent  and  with  a length  of  about  17  meters,  the  ground  wires 
actually  decrease  the  low  angle  radiation  by  perhaps  7 or  8 percent. 

The  next  Figure,  9,  gives  the  corresponding  data  for  a lower  frequency, 
namely 3^Hz,  with  the  same  three  values  of  earth  conductivity.  These  data 
show  that  a)  at  the  lower  frequency,  the  possible  enhancement  is  less, 
being  perhaps  3?  for  10  milllmho/m  earth,  6%  for  5 mi  11  imho/m  earth  and 
perhaps  18%  for  1 mi  11  imho/m  earth,  b)  the  ground  screen  length  to  accomplish 
the  maximum  ehancement  Is  longer,  in  the  vicinity  of  19  to  20  meters  with  the 
lowest  va i ue  of  earth  conductivity. 
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FIELD  ENHANCEMENT  AS  A FUNCTION  OF 
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The  next  Figure,  10,  gives  the  corresponding  result  for  30  MHz,  and 

of  course  shows  both  higher  possible  enhancements  and  more  drastic  variations 

In  the  enhancement  with  length.  Again,  the  variation  in  the  enhancement  with 

length  is  less  with  the  higher  conductivity  values.  And  the  data  shows  that 

-2 

enhancements  of  about  9 to  19%  are  possible  with  conductivity  10  mhos/meter. 

The  foregoing  data  were  obtained  on  the  assumption  and  on  imposition  of 
a traveling  wave  distribution  of  current  on  the  ground  wires,  largely  having 
a propagation  constant  equal  to  Coleman's  approximation  p.  10,  eqn.  6. 

Coleman's  approximation  is  valid  only  for  the  wi res  essent ia 1 1 y lying  in  the 
interface.  If  the  wires  have  significant  insulation  around  them  or  if  they 
are  somewhat  above  the  interface,  the  propagation  constant  changes  significantly. 
To  document  the  effect  of  such  differences  in  propagation  constant,  the  calcula- 
tions of  the  type  given  in  Fig.  10  (30  MHz)  were  repeated  with  different  (Fig.  11) 
assumed  values  for  the  propagation  constant  of  the  traveling  wave  on  the 
ground  wire.  A number  of  different  propagation  constants  were  assumed,  to 
show  the  effect  of  both  a change  in  the  real  part  (attenuation)  and  the 
imaginary  part  (phase  constant)  of  the  propagation  constant.  Of  course,  with 
the  higher  attenuation  constants,  the  enhancement  tends  to  reach  a maximum 
at  some  length  and  to  level  off  thereafter.  For  the  phase  constants  smaller 
than  that  of  the  Coleman  formula,  the  length  that  shows  the  maximum  enhancement 
tends  to  be  longer  and  if  the  attenuation  is  small,  the  enhancement  levels 
larger. 

This  dependence  of  the  results  on  the  actual  wire  propagation  constants 
Is  important  from  a practical  viewpoint  since  the  positioning  and  Insulation 
thickness  of  the  ground  wires  is,  at  least  to  some  extent,  under  the  control 
of  the  system  designer. 





FIG.  11 

FIELD  ENHANCEMENT  AS  A 
FUNCTION  OF  LENGTH  AT  30  MHZ 
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Of  course,  when  the  WF-LLL2AP  program  Is  used  In  the  calculation,  no 
assumption  regarding  propagation  constant  or  current  distribution  on  the  wires 
Is  necessary.  The  program  output  gives  the  current  values  at  the  center  of 
each  segment  length  used  In  the  calculation.  From  these  data  It  is  possible 
to  obtain  an  approximation  for  the  propagation  constant  on  the  ground  wires. 

This  was  done,  at  a frequency  of  10  MHz  for  two  values  of  conductivity 
*3  “2 

(10  3 and  10  mhos/m)  wi  th  a wi  re  of  3 . 77mm  rad  ius  located  2 . 1 rad  i i above  the  i nter- 
face.  In  this  way  it  was  found  that  the  Coleman  result  gives  a propagation 
constant  that  is  reasonably  close  to  the  apparent  propagation  constant  as 
inferred  from  the  calculated  current  distribution. 

Another  important  variable,  to  some  extent  under  the  control  of  the 
designer,  is  the  type  of  termination  at  the  end(s)  of  the  ground  wire.  The 
easiest  alternative  in  practice  is  of  course  to  leave  the  far  end  open 
circuited.  But  for  low  conductivity  earth  and  relatively  short  ground  wires, 
this  may  not  be  a desirable  alternative.  The  reason  is  that  the  ground  wire 
system  may  tend  to  exhibit  resonant- impedance  conditions  for  certain  lengths. 

The  system  effects  of  such  conditions  Is  suggested  by  the  data  in  Fig.  12; 
which  were  obtained  by  running  the  computer  program,  modified  WF-LLL2AP. 

Figure  13  shows  similar  type  effects  when  the  approximate  calculation  was 
carried  out  on  the  assumption  of  an  open  circuit  termination  on  the  ground 
wire.  The  effect  shows  up  in  a striking  way  because  of  the  internally 
imposed  condition  that  the  total  current  at  the  input  of  the  ground  wires  is 
equal  to  the  total  current  at  the  base  of  the  vertical  antennas.  The  main 
thing  this  data  shows  Is  that  ther**  can  be  resonant  effects  on  the  ground 
wires  under  the  conditions  (low  earth  conductivity)  in  which  the  ground  wires 
are  the  most  important.  The  input  impedance  indicated  in  the  WF-LLL2AP  program 
also  varies  wildly  as  a function  of  ground  wire  length. 
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FIG.  13 

FIELD  ENHANCEMENT  AS  A FUNCTION  OF  LENGTH 
AT  10  MHZ,  APPROXIMATE  CALCULATION. 
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Short  circuiting  (low  impednnce)  terminations  on  the.  ground  wires  can 
give  similar  effects,  although  In  practice,  such  low  Impedance  terminations 
would  be  difficult  to  achieve  In  practice  if  the  earth  conductivity  is  small. 

Conclusions 

After  some  months  of  effort,  we  have  concluded  that  it  is  possible  to 
employ  a method  of  moments  program  incorporating  a Sommerfeld  type  of  dipole 
field  solution  to  evaluate  the  effects  of  wire  ground  systems  in  the  HF  band. 
However,  with  a computer  of  the  type  available  at  Purdue  University  (CDC  6500) 
it  is  so  far  only  marginally  so.  It  is  marginal  1)  because  the  evaluation  of 
the  "self-impedance"  and  nearby  "mutual  impedance"  type  terms  requires  so 
much  time  and,  2)  because  of  the  limitations  on  the  system  size  imposed  by 
the  difficulties  and  the  time  required  to  invert  the  associated  large  order 
system  impedance  matrix. 

By  a comparison  of  the  results  of  the  digital  computer  output  with  the 
approximate  calculation,  we  conclude  that  the  assumption  of  a transmission 
line  like  current  distribution  having  a propagation  constant  calculated  in 
standard  fashion  (such  as  Coleman  result)  can  be  employed  to  give  useful 
information  on  the  radiation  effects  of  wire  ground  systems. 

In  general,  the  ground  wires  may  enhance  the  radiation  at  all  forward 
and  upward  angles.  However,  the  amount  (and  even  whether  the  field  is  increased 
or  decreased)  of  the  enhancement  depends  on  the  ground  wire  length  and  propagation 
constant.  For  low  conduct  ivi  ty  earth,  resonance  effects  occur  with  certa  in  lengths  . 

For  low  conductivity  earth,  enhancements  of  about  20%  appear  to  be 
possible,  but  the  results  are  length  and  frequency  sensitive. 

Wire  lengths  that  are  long  enough  or  properly  terminated  so  as  to  give  a 
traveling  wave  distribution  on  the  ground  wires  have  potential  bandwidth 
advantages . 
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From  a systems  point  of  view,  the  fact  that  the  designer  has  some  control 
over  the  ground  wire  propagation  constant  (through  insulation  thickness)  and 
over  the  detailed  type  of  grounding  connections  may  be  of  vital  importance  to 
the  enhancement  of  low  angle  radiation. 

On  the  negative  or  discouraging  side,  it  must  be  noted  that  in  no  case 
studied  so  far  can  the  field  enhancement  be  said  to  really  alter  the  basic  "null 
at  the  horizon"  effect;  for  even  an  increase  by  a factor  of  two,  which  might  be 
possible,  still  leaves  the  field  near  the  horizon  at  quite  a small  value.  Moreover, 
the  patterns  show  an  increase  of  field  strength  in  the  whole  forward  quadrant 
and  it  is  not  yet  clear  that  the  power  represented  by  these  field  increases 
comes  from  power  that  might  otherwise  have  been  dissipated  in  the  earth  or 


whether  the  additional  power  would  have  to  be  supplied  by  the  transmitter. 
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METRIC  SYSTEM 


BASE  UNITS: 

Quantity 


Unit 


length 

mass 

time 

electric  current 
thermodynamic  temperature 
amount  of  substance 
luminous  intensity 

SUPPLEMENTARY  UNITS: 

plane  angle 
solid  angle 

DERIVED  UNITS: 

Acceleration 

activity  (of  a radioactive  source) 

angular  acceleration 

angular  velocity 

area 

density 

electric  capacitance 

electrical  conductance 

electric  field  strength 

electric  inductance 

electric  potential  difference 

electric  resistance 

electromotive  force 

energy 

entropy 

force 

frequency 

illuminance 

luminance 

luminous  flux 

magnetic  field  strength 

magnetic  flux 

magnetic  flux  density 

magnetomotive  force 

power 

pressure 

quantity  of  ele<  trie  ity 
quantity  of  heat 
radiant  intensity 
specific  heat 
stress 

thermal  conductivity 
velocit) 

viscosity,  dynamic 

viscosity,  kinematic 

voltage 

volume 

wavenumber 

work 


metre 

kilogram 

second 

ampere 

kelvin 

mole 

candela 


radian 

steradian 


metre  per  second  squared 

disintegration  per  second 

radian  per  second  squared 

radian  per  second 

square  metre 

kilogram  per  cubic  metre 

farad 

siemens 

volt  per  metre 

henry 

volt 

ohm 

volt 

joule 

joule  per  kelvin 

newton 

hertz 

lux 

candela  per  square  metre 
lumen 

ampere  per  metre 

weber 

tesla 

ampere 

watt 

pascal 

coulomb 

loule 

watt  per  steradian 

loule  per  kilogram  kelvin 

pascal 

watt  per  metre-kelvin 

metre  per  second 

pascal-second 

square  metre  per  second 

volt 

cubic  metre 
reciprocal  metre 
joule 


SI  PREFIXES 


Multiplication  Factors 
1 000  000  000  OCX)  10" 

1 000  ooo  ooo  - 10” 

1 ooo  ooo  - 1 0* 

1 000  - 10’ 
ioo  - inJ 
10  = 10’ 

0 1 10-' 
001  10-' 

0 001  10-’ 
0 1100  001  10  * 

o OCX)  oiH)  ooi  io  • 

o ooo  non  non  ooi  io  " 

0 000  000  OOO(HK)  001  10  " 

o ooo  ixx)  ooo  txio  txio  ooi  to  •* 

• Tn  hr  avoided  where  possible 


SI  Symbol  Formula^ 


m 

kg 

s 

A 

K 

mol 

cd 


rad 

sr 


mis 

(disintegration  )/s 

rad/s 

radis 

F 

m 

kg/m 

A-s/V 

S 

A/V 

H 

V m 
V-s/A 

V 

W/A 

V 

V/A 

W/A 

) 

N-m 

)(K 

N 

kg-m/s 

Hz 

(cycle)/s 

lx 

im/m 

Im 

cd/m 

cd-sr 

Wb 

A/m 

V-s 

T 

Wb/m 

A 

W 

V* 

Pa 

N/m 

C 

) 

A-s 

N-m 

Pa 

W/sr 

1/kg-K 

N/m 

V 

W/m-K 

m/s 

Pas 

m/s 

W/A 

m 

(wave)/m 

N-m 

Prefix 

.SI  Sym 

tera 

T 

gig* 

megs 

(f 

M 

k 

kilo 

hecto* 

h 

deka* 

da 

deci* 

d 

centl* 

c 

milli 

m 

micro 

M 

nano 

n 

pico 

femto 

r 

si  to 

a 

